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Abstract. In this paper we calculate the emission of gravity waves by the binary 
pulsar in the framework of five dimensional spacetime. We consider only spacetimes 
with one compact extra-dimension. We show that the presence of additional degrees of 
freedom, especially the 'gravi-scalar' leads to a modification of Einstein's quadrupole 
formula. We compute the induced change for the binary pulsar PSR 1913+16 in 
the simple example of a 5d Minkowski background. In the example of a cylindrical 
braneworld it amounts to about 20% which is by far excluded by present experimental 
data. 
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1. Introduction 

During the last couple of years, the possibility that our 3+1 dimensional Universe might 
be a hypersurface in a higher dimensional spacetime has been investigated. This can 
be consistent with relatively 'large' extra-dimension if only gravitational interactions 
are propagating in the full spacetime, while gauge interactions are confined to the 3 + 1 
dimensional 'brane' PQEj- The motivation for these 'braneworlds' is twofold. First of all, 
super string theory, which is consistently formulated only in 10 dimensions, predicts such 
branes as endpoints of strings, 'D-branes' [3]. Secondly, allowing for n extra dimensions 
of size L, the observed, four dimensional Planck scale, M4, relates to the fundamental 
Planck scale of the underlying string theory, Mp via M 4 = Mp(LMp) n / 2 . Hence for a 
sufficiently large length L a small fundamental Planck scale, Mp ~ TeV can lead to the 
observed effective Planck scale, M 4 = I/^AttG^ ~ 3.4 x 10 18 GeV, where G 4 denotes 
Newton's constant PQ E]- A short calculation shows that for n > 2 the required size 
of the extra dimensions is less than about 1mm at which scale Newton's law has not 
been extensively tested so far. Here and in the following we use natural units to that 
h = c=l. 

Motivated by this new attempt to solve the hierarchy problem, people have started 
to study the modified gravitational theory obtained on the brane, when starting from 
Einstein's equation in the higher dimensional 'bulk'. Many astrophysical and laboratory 
consequences of this idea have been investigated and found to be consistent with present 
bounds [1] . For simplicity, most of the work has concentrated on 1 single extra dimension 
for simplicity. We also do so in this work. In Refs. 010] the modifications of cosmo logical 
solutions are derived. Following that, a lot of work has been devoted to cosmo logical 
perturbation theory in the context of brane worlds (see Refs. [7j to j2H| and references 
therein). To do this, the authors have simply linearized the 5-dimensional Einstein 
equations around a cosmo logical background solution. In this paper we show that, at 
least for compact braneworlds, which allow a zero-mode like it is known from Kaluza — 
Klein theories, the five dimensional linearly perturbed Einstein equations are in conflict 
with observations. 

For this we will make use of the following very generic issue: In d + 1 dimensions, 
the little group of the momentum state (p, 0, • • ■ , 0,p) of a massless particle is SO(d— 1) 
(see e.g. [2I])- For d = 3, the irreducible representations of 5*0(2) together with parity 
lead to the usual two helicity states of massless particles, independent of their spin. In 
d = 4 dimensions however, a massless particle with spin 2 like the graviton has the 5 spin 
states coming from the five dimensional tensor representation of SO (3). Projected onto 
a 3-dimensional hypersurface, these five states become the usual graviton, a massless 
particle with spin 1 usually termed 'gravi-photon' and a massless particle with spin 0, 
the 'gravi-scalar'. These particles couple to the energy momentum tensor and should 
therefore be emitted by a time dependent mass/energy configuration like a binary system 
of heavenly bodies. This problem has not been discussed in the literature before. In 
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analogy to the four dimensional quadrupole formula for graviton emission, 

q 

~ ^ = 45^ Q ^iQw (!) 
one might naively expect a formula for dipole emission of for the gravi-photon, like 

- E oc -SiSi and (2) 
c 3 

—E oc —MM (3) 

c 

for the monopole emission due to the gravi-scalar respectively. Here S is the center of 
inertia, 



Sj = J d?x p(x)xj 



so that S = P and M is the mass of the system. Contrary to Q, the momentum, S = P 
and the mass are conserved to lowest order in G. Therefore, the contributions (J2J) and 
(J3J) vanish and the quadrupole formula, which is so well tested with the binary pulsar 
PSR 1913+16 that Hulse and Taylor have obtained the Nobel price of 1993, is 
maintained to sufficient accuracy. 

In this work we show that Eq. (J2J) is indeed verified and there is no emission of gravi- 
photons to lowest order in G. However, Eq. Q is not. In formulating it we omitted 
another scalar quantity, the trace of the second moment on the mass distribution, 
I = J px^x^d 3 x. As we shall derive here, this term leads to an emission of gravi-scalars 
which induce an energy loss of the same order as the usual quadrupole term. We will 
see that for the binary pulsar PSR 1913+16 this modifies the prediction from Einstein's 
quadrupole formula by about 20% in blatant contradiction with the observations which 
confirm the formula within an error of about 0.5%. 

In the next section we derive the result announced above. In Section III we discuss 
its consequences for models with one extra-dimension. 

Throughout we use capital Latin letters for five dimensional indices 0,1,2,3,4, 
Greek letters for four dimensional ones and lower case Latin letters for three dimensional 
spatial indices. The spatial Laplacian is denoted by, S^didj = A, and spatial vectors 
are indicated in bold face. 



2. The modified quadrupole formula for spacetimes with one 
extra— dimension 

In order to be specific and for simplicity, we consider a 4 + 1 dimensional cylinder, with 
a rolled up fourth spatial dimension of length L. Since this length which is of the order 
of micrometers or smaller, is much smaller that the scale of the 3 + 1 dimensional system 
which is emitting gravity waves, we shall only consider states which are zero-modes with 
respect to the fourth spatial dimension. 
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We choose a gauge such that the perturbed five dimensional line element is given 
by [H] 

ds 2 = - (1 + 2m)dt 2 - 2Y, i dtdx i - 2Bdtdy + 2S i dx i dy 
+ [(1 - 2$)<% + 2Hij] dx i dx j + (1 + 2C)dy 2 
= (vab + h AB )dx A dx B . (4) 

Here Ej and Si are divergence free vectors, iiy is a traceless, divergence free tensor 
and tjab is the five dimensional Minkowski metric. In the appendix we show, that the 
variables defined above are gauge-invariant combinations of the most general metric 
perturbations. There we also see that the gauge choice made in Eq. (J1J) fixes the gauge 
completely. The extra-dimension is parameterized by the co-ordinate y, < y < L for 
our cylindric spacetime. We are working within linear perturbation theory and therefore 
require \Hab\ ^ 1- The source free linearized Einstein equations for this geometry can 
be reduced to three wave equations, 

(5) 
(6) 
(7) 

The other variables are then determined by constraint equations, 

AC = (2A + 3d 2 )<5> (8) 
# = $ - C (9) 
AB = - 6d t d y $ (10) 
d v Ei = - <9 t £ 4 (11) 

(for more details see the Appendix and |26j). The important point is that the above 
wave equations describe simply the five propagating modes which we also expect from 
our group theoretical argument above. We now consider a four dimensional, Newtonian 
source given by 

T M = T 44 = (12) 
Too =p(x,t)6(y) (13) 
T 0l = (diV + Vi)(x,t)5(y) (14) 







- A)$ 


= 






- A)£j 


= 


(d t 2 




- A)H i3 


= 



T ■ 



+ i(^ + W + n, 



(x,t)6(y) (15) 



where V{ and Ilj are divergence free vector fields and ITj is a traceless divergence free 
tensor in three spatial dimensions. Eqs. (|12j) to (|15p represent the most general energy 
momentum tensor on the brane with vanishing energy flux into the bulk. However, being 
interested in the binary pulsar, we shall restrict our attention to Newtonian sources, i.e. 
sources with p = \T 00 \ 3> |T 0i | ~ pV , |T 00 | ^> \T^ ~ pV 2 , where V 1 is a typical 
velocity of the system. 
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Since we are only interested in the zero-mode with respect to the fourth spatial 
dimension, we may integrate over the fifth dimension. The four dimensional linearized 
Einstein equations with source term then become (see Appendix) 



scalar perturbations: 

A(2$ -C) = Sudp (16) 

d t C - 2d t <$> = 8ttG 4 v (17) 

d t 2 {-C + 2<t>) = 8ttG 4 (p + ^AIl) (18) 

$ - - C = 8nGJl (19) 

AB = , d t B = (20) 

-A(2$ - + 3<9 t 2 $ = (21) 
vector perturbations: 

^ASi = torG&i (22) 

d t Ei = torGJU (23) 

A£ 4 - d t 2 £ = (24) 
tensor perturbations: 

-AH i:i + d t 2 H i:j = 87rG 4 n ii . (25) 



where G 4 = (ir/2)L 1 G§ is the four dimensional Newton constant, M 4 = (2ir 2 G5/ 1 L) 1 / 2 . 
Here we have used the relation G d l = M$~ 2 vo\(S d - 2 ) = Mf~ 2 2ir i r /T(^) between 
the (i-dimensional Newton constant and Planck mass. This relation ensures that Gauss' 
law is true in any dimension. Einstein's equations in an arbitrary number of dimensions 
are Gab = 2M 2 L ~ d TAB- Note that the equations for the zero- mode are very similar 
to the Kaluza-Klein approach. Both B and Si are completely decoupled. In Kaluza- 
Klein theories they play the role of the electromagnetic field which couples to the current 
J M oc T^ 4 which we have set to zero in this work. In a coordinate system where the brane 
is fixed (no brane bending) the only difference for the zero-mode, in our 'test-brane' 
analysis is the relation between the four- and five-dimensional gravitational constant. 
From Eqs. (|16|) to (|25|) we can derive three wave equations with source term: 



(<9 2 -A)$ = -87rG 4 (p-An) (26) 

(<9 2 -A)£ 4 = 87rG 4 (IL, + 2 Vl ) (27) 

d t Ei = 87rG 4 IL, (28) 

(d 2 - A)H l3 = 87rG 4 U i:j , (29) 



The other variables are again determined by constraint equations. 

Eq. ()29|) leads over well-known steps to the quadrupole formula for Newtonian 
sources. The retarded solutions of the wave Eqs. ()26|) and (|2~7jl determine Ej and $ 
far away from the source. It is easy to see that the a priori dominant terms p in (|26J) 
and Vi in (12 7|) lead exactly to the terms predicted in Eqs. © and (J2J). Hence their 
contribution to the energy emission vanishes (to lowest order) by mass and momentum 
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conservation. We therefore concentrate on the terms coming from the stress tensor. 



Energy momentum conservation (T' 



0) gives 



A(AII + ^P) = ^didjT ij = ^<9 2 p 



so that All = l(U - P), where 



U(t,x) 



p(t,x') 



d 3 x' . 



(30) 



(31) 



Air J5R3 |x — x', 

Taking the divergence of (|15|) and using the above expression for All we find 

AIL j = 2d i (Tj-S i j U) , (32) 

so that, to order 1/r, each component of the vector potential Ilj is the integral over 
a divergence which vanishes for a confined source. Hence Ilj oc r -2 ; together with 
()28|) this implies that also is of order r -2 . The energy flux from Sj is expected 
to be proportional to t Q j oc (9jEj)(9[ i E :7 -]) oc 1/r 4 (the bracket [..] indicates anti- 
symmetrization) which shows that we have no energy emission from the gravi-photon 
far away from the source. 

The situation is different for the gravi-scalar, $. There the energy flux is 
proportional to £o« dt§di§. To order 1/r = l/|x|, the scalar field $ is given by 



*{t, x) = 9± [ \3(U(t', x') - P(f , x')) - 2p(t', x') 



d 3 x'. 



(33) 



Here f = t — |x — x'| is retarded time. The last term does not contribute to t 0i since 
dt J p = M = 0. Furthermore, p = <9jT° l so that 



4nU{t, x) 



djT^jt,^) 
x — x'| 

Oi 



d 3 x' 



T U4 (t,x')a;- 



rttV 



0; 



T 0l (t,x') 



(34) 
(35) 



Here <9 4 ' denotes differentiation w.r.t. x' while di is differentiation w.r.t. x. In other 
words, U is a divergence and therefore does not contribute in the integral ("3~3~j) . 
Finally, we use the identity (see [2*7] ) 

^//^ 

which is a simple consequence of energy momentum conservation and Gauss' law. All 
this leads to 



d' J T, 



$(t, x) = ^ / p (t', x')x' W + 0{Aj 



and 



2r 
"27 



(36) 
(37) 



Here n* 



x. l /r. The energy emission from the gravi-scalar is then given by 



'dE\ 



scalar 



lim / r dttn l t 0i (m) 
r ^°° Js 2 



oc lim 

r— >oo 



eft 3 



px 2 d 3 x 



(38) 
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For the last equality sign we have used the results (JHtj|) and (|37|) and integrated over the 
sphere. 

To obtain the proportionality factor which is still missing in Eq. (j38|) . we write the 
five dimensional action up to second order in the scalar metric perturbations 

s = ^cjl dy L dAx ^ R 

= TfiV L*x(l + hi)(RM + RV) 

107rG4 jk 4 2 

= J d^x-d^ipd^ip + a total derivative . (39) 

Here h\ = 2(\E' — 3$ + C) is the trace of the metric perturbations and and R^ 
are the first and second order (scalar) perturbations of the five dimensional Riemann 
scalar. The last equation is just the ansatz for a usual, canonically normalized scalar 
field action. The source free Einstein equations for the zero-mode can be obtained by 
varying this action. Therefore, the energy flux of the mode is given by toi — doypdiip. A 
straight forward but somewhat tedious calculation of R^ and R^ 2 \ in which we use the 
homogeneous constraint equations to express all scalar perturbations variables in terms 
of $, gives finally 

/ 21 21 • 

(p — \ — — $ so that t 0i = — — $<9i$ . 

Inserting this proportionality factor in Eq. (J3%|) leads to 

(40) 



(dE\ 


21G 4 


~d 3 [ 2 


\dt) 


scalar « 


Tw3 / , P x 

dt 6 



A textbook calculation for a binary system in Keplerian orbit, like it is presented 
e.g. in Ref. j2Zj, now gives for the energy loss averaged over one period of the system 

_/dE\ 21GjM^Ml(M 1 + M 2 )g(e) 
\dt/ ~~ 4a 5 (l-e 2 ) 7 /2 [ > 

with g(e) = e 2 (l + e 2 ) . (42) 

Here M\p are the two masses of the system, a is the semi-mayor axis and e is the 
eccentricity of the orbit. Comparing this with the quadrupole formula for a Keplerian 
binary system which can be found e.g. in [2Z|), we see that (|41j) agrees with it upon 
replacing ^1 by ^HiM w i t h f = 1 + (73/24)e 2 + (37/96)e 4 . Inserting the value 
e = 0.617 for PSR 1913+16 given in one finds that this new contribution amounts 
to 19.9% of the ordinary quadrupole prediction. Therefore, the slowdown of the period 
of PSR 1913+16 should be 

T|tens+ S cai = - 2.88 x 10~ 12 instead of (43) 
T| tens = - 2.4024 x 10- 12 . (44) 

But the experimental value is in perfect agreement with the quadrupole formula (|44jl 
and contradicts such a huge correction: 

f\ cxp = -(2.408 ± 0.01) x 10~ 12 . (45) 
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Our result is not entirely self-consistent in that we assume that the masses Mi 
and M 2 inferred for the two neutron stars are the same as in ordinary four dimensional 
gravity. This is not necessarily the case and probably also depends on the method 
by which they are determined. In principle we would have to specify two additional 
measurements and check whether different masses M[ and M' 2 can be found so that 
all the experiments agree with the results from four dimensional gravity for the masses 
Mi and M 2 . This procedure has been undertaken in the work on general scalar-tensor 
theories of gravity in Ref . |2H] . Clearly, it would be a very surprising coincidence if such 
masses could be found. 

Finally, we want to mention that other tests, like for example light deflection around 
the sun most probably also lead to deviations from the four dimensional result. 



3. Conclusion 

Our calculations show clearly that a 4 + 1 dimensional cylindrical spacetime which 
satisfies the five dimensional Einstein equations contradicts observations. This finding 
is independent of the size of the extra dimension. Note that we did only take into 
account the Kaluza-Klein zero-mode and hence the size L of the extra dimension just 
enters in the relation between the five and four dimensional gravitational constants, 
G4 = G 5 7r/2L. It is important to note at this point that the zero- mode in general of 
course depends on the background geometry and thus on the warp factor (it is a solution 
to the five dimensional d'Alambertian, D 5 with respect to the background geometry). 
In our case, the background geometry is trivial. If the warp factor depends on the size 
of the extra-dimension this will of course also imprint itself on the zero-mode. 

In our calculation we have considered a simple cylindrical spacetime. But our result 
would not alter substantially if we would consider a more realistic, warped spacetime 
which corresponds to an expanding universe. The mayor exception are warped models 
with non-compact extra-dimensions like the Randall Sundrum II model [2H], where the 
scalar gravity mode is not normalizable and therefore our discussion does not apply. But 
in all models with compact extra- dimensions, for example in all models with two branes, 
like Randall Sundrum I [3U] or the ekpyrotic universe the mode discussed here is 
normalizable. Its coupling is determined by the five dimensional Einstein equations. 
We therefore expect that (besides all the other problems of scalar-tensor theories of 
gravity 28, ) this mode be excited by a system like the binary pulsar and lead to 
deviations from Einstein's quadrupole formula. Even though the exact amount of change 
will depend on the details of the model, generically we can expect it to be of the same 
order of magnitude. The expansion of the Universe, which determines the jump of 
the extrinsic curvature between the two sides of the brane is not relevant on scales of 
a binary pulsar and so we expect a realistic braneworld model to lead to very similar 
results. The possibly more complicated geometry of the extra-dimension may somewhat 
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affect the numerical result but as long as the zero-mode discussed here exists (i.e. as 
long as it is normalizable) the problem will remain. 

In order to find consistent four dimensional Einstein gravity on the brane, the five 
dimensional Einstein equations have to be modified somehow at very low energy. One 
possibility might be that the gravi-scalar (p acquires a small mass which is larger than 
the frequency of the binary pulsar, > 1/T ~ 3.6 x 10 _5 Hz ~ 2.4 x 10~ 20 eV (e.g. via 
the Goldberger-Wise mechanism [32 ). Such a tiny mass would already suffice to inhibit 
the field p to be emitted by the binary pulsar PSR 1916+13. To interpret the gravi- 
scalar physically it is useful to notice that in source free space, for the zero modes we 
have 2$ = C. But C determines the size of the extra dimension, L + 5L = / L (1 + C)dy. 
Therefore, if one wants to stabilize the size of the extra dimensions one needs a positively 
curved potential, hence a mass, for C and thus for $. This argument is well known in 
string inspired Kaluza-Klein theories. There, the size of the extra dimensions is given by 
the string scale, i and T-duality, the duality of string theory under the transformation 
r — > £ 2 /r, indicates that a potential for 5L and hence for C should have a minimum at 
L = £, C = 0, in other words a positive mass term, see e.g. Some non-gravitational 

(stringy) correction must provide this mass. 

But even if we set C = 0, the scalar perturbation equations (fTB^) to (}2Tj) do 
not reduce to the usual four dimensional perturbation equations: using (fTSj) and the 
Laplacian of (fl"9|) to eliminate <9 f 2 $ and A($ — \1>) in (|2"Tf . we obtain 

A$ = 8ttG 4 ^P (46) 

Comparing this with Eq. (|16|) we find that these equations are only compatible when 
p = 3P, in other words, if the matter is relativistic. This shows that the problem 
discussed in this paper goes beyond the radion problem: even if the radion C is fixed, 
we do not recover the four dimensional Einstein equations. 

Moreover, if non-gravitational corrections influence the evolution equation of the 
gravi-scalar, they may very well also introduce modifications to the interactions on the 
brane which are usually inferred from the five dimensional Einstein equations, like the 
Israel junction conditions [SI], or certainly the gravitational perturbation equations as 
derived in Refs. [Zj to 

As mentioned above, the situation is mitigated if the geometry is warped, like in 
the Randall- Sundrum model (2H!- In such models, the fifth dimension may be non- 
compact and the gravi-scalar zero-mode can be non-normalizable in the sense that the 
integral / \J~\g\ip 2 dy = oo. Then the gravi-scalar cannot be excited as it has infinite 
energy. Within the language of classical relativity, we may simply say that the 'gravi- 
scalar mode' turns out to be large for some values of y and linear perturbation theory 
is therefore not applicable. 

Let us finally summarize our main conclusion: If a model with extra dimensions 
is such that the gravi-scalar coming from Einstein's equations in five dimensions is 
normalizable, which is always the case if the extra-dimensions are compact, non- 
gravitational interactions need to intervene in order to make the model compatible 
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with observations. They have to give the gravi-scalar a mass. In such models 
therefore, we cannot simply use the five dimensional Einstein equations to infer the 
four dimensional equations of motion. Especially, we cannot apply five dimensional 
gravitational perturbation theory to, e.g., a cosmological solution in order to learn 
about four dimensional cosmological structure formation. The perturbation theory 
derived in Refs. [7j to [2H] can at best be applied to models where the gravi-scalar 
is not normalizable. 

Clearly, the situation does not improve if there are more than one extra dimensions. 
In this case there are several gravi-scalars which may all contribute a similar amount to 
the emission of gravity waves. 
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Appendix A. 

Gauge— invariant perturbation equations for a 5d Minkowski braneworld 

We consider a five dimensional Minkowski bulk with background metric (tjab) = 
diag(— 1, 1, 1, 1, 1, ). The most general perturbed metric is then of the form 

ds 2 = - (1 + 2A)dt 2 - 2S l dtdx i + (<J y + h i3 )dx l dx^ - 2Bdydt 



The spatial vectors Si and Ei can be decomposed into scalar (spin zero) and vector 
(spin one ) components and the spatial tensor hij can be decomposed into scalar, vector 
and tensor (spin two ) components; 



Here E i} Fi and Si are divergence free vectors and is a divergence free, traceless 
symmetric tensor. To first order, scalar vector and tensor perturbations evolve 
independently. Spatial indices can be raised and lowered with the background metric 



Let us now study the transformation of the defined perturbation variables under 
infinitesimal coordinate transformations (gauge transformations) in the bulk. We 
consider an infinitesimal coordinate transformation 





Si = ViS + Si, Ei = ViE + Ei, 

= 2i?7jj + 2Fij, Fij = V(ji^) + H^, Fi = VjF + Fj. 



(A.2) 
(A.3) 




x 



A 




(A.4) 



where we set 




A 



i 




(A.5) 
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The three scalar (T, L, L 4 ) and one vector type (L ) gauge-transformation allow us to 
gauge to zero three scalar and one vector variable. Instead of choosing a specific gauge 
(like in the main text of this paper) we shall define gauge invariant combinations and 
express the first order Einstein equations in terms of these. This is always possible as a 
consequence of the Steward- Walker Lemma |35j . 

Under the above coordinate change, the geometrical perturbations transform 



according to the Lie derivative of the background metric, qab 9ab + L^ab 

A^A + dtT, Si^Si- d t L h (A.6) 

H^H, F^Fi + Li, (A.7) 

Ha -+H ijt B^B- d t U + d y T, (A.8) 

E l ^E i + d y L h C^C + d y L 4 , (A.9) 

3 -»• S + T-d t L, F^F + L, (A.10) 

E —> E -{- d y L + L A . (A.ll) 

We can therefore define the following four scalar and two vector perturbation variables 
which are gauge invariant, 

* =A-dt(S + dtF), (A.12) 

$ = -H, (A.13) 

B =B-dy(S + d t F) + d t (E — d v F), (A.14) 

C =C-d y (E-dyF), (A.15) 

E i = S i + dtF i , (A. 16) 

£ i =E i -d y F i . (A. 17) 



The tensor variable Hij is gauge invariant since there are no tensor type gauge 
transformations. 

The stress energy tensor defined in Eqs. (fT2|) to (fTK|) is a first oder perturbation 
with vanishing background component and it is therefore gauge invariant by itself. This 
is a specialty of the Minkowski background. The much more complicated general gauge- 
invariant variables for a generic evolving bulk can be found in Ref. 

From Eqs. (jA.7|) and (jA.ll|) it is clear that the gauge is completely fixed by setting 
F = S = E = and Fj = 0. This 'generalized longitudinal gauge' is precisely the gauge 
choice adopted in the main text of this paper. 

We can now go on and compute the components of the curvature in terms of our 
perturbation variables. In Ref. [2H] this is done for the generic evolving case. Here we 
repeat for convenience the expressions for the Christoffel symbols and for the Einstein 
tensor in our simple case. In contrast to the generic situation, the Einstein tensor is 
gauge invariant in for a Minkowski bulk. 

Christoffel symbols 

T Q ° = d t A, T°=diA, (A. 18) 
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r ij - V(iSj) + -d t hij, 


1 04 




(A.19) 


r,° 4 = \ (d y S t + d t B + d t Ei) , 


1 44 


= 9,5 + d t C, 


(A.20) 


= &A - d t s\ 


r * 

1 Oj 


= \ {#Si - djS* + d t h)) , 


(A.21) 


r j i k = l(d j hi + d k h)-#h jk ), 


1 04 


2 v " / 


(A.22) 


rt — - (d h?- - d l E + d-E') 

j4 2 V V 3 3 ] ) 


i 1 44 


= d y E i - d l C, 


(A.23) 


r 4 = d y A - ^B, 


r 4 

1 Oi 


1 , 

= - (d y Si - d t B + 8^) , 


(A.24) 


T tj = ~2 9 y hi 3 + 
= diC, 


r 4 

1 04 

r 4 

1 44 


= d t C, 

= 0yC, 


(A.25) 
(A.26) 



Here X(jXj) = ^(XiYj + XjYi) indicates symmetrization in the indices i and j. 
Einstein tensor 

Over the Ricci tensor one can now compute also the Einstein tensor. This is a somewhat 
lengthy but straight forward calculation. Also here, the more complicated results for an 
evolving bulk can be found in Ref. 



Goi 
n , . 



G 



0-1 



a 



A(-C + 2$) + 3d 2 $, 



di 



d t (2$-C)--d y B 



^dy {dyZi + d^} + ^ASi 



[ASy - didj) (* - $ + C) + <y<9 t 2 (-C + 2$) + cP(* - 2$) 

- dtdyB] + dtdliEfl + dyd^ + {d^ - A - fyHij, 



AB + 3dtd y $, 
d y (2<$> - *) + + ^t^S, + «9 2 £ - A£ 



G 44 = A(tf-2$) + 3^ 2 $ . 



(A.27) 
(A.28) 

(A.29) 
(A.30) 

(A.31) 
(A.32) 



The source free Einstein equations (J3J) to are now simply obtained by setting 
Gab = 0. Integrating the five dimensional Einstein equation, Gab = 2M% 3 Tab over 
the extra dimension y, with the ansatz given in Eqs. (|12j) to (|15p for the stress energy 
tensor one obtains the 4-dimensional Einstein equations given in Eqs. (jl6j) to (|25jl. 



References 



[1] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B429, 263 (1998). 

[2] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B436, 257 (1998). 

[3] J. Polchinski, String Theory, Vol I & II, Cambridge University Press (1998). 



Testing extra dimensions with the binary pulsar 



13 



[4] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Rev. D59, 086004 (1999). 

[5] P. Binetruy, C. Deffayet, and D. Langlois, Nucl. Phys. B565, 269 (2000). 

[6] P. Binetruy C. Deffayet, U. Ellwanger, and D. Langlois, Phys. Lett. B 477, 285 (1999). 

[7] C. van de Bruck, M. Dorca, R. Brandenberger, and A. Lukas, Phys. Rev. D 62, 123515 (2000). 

[8] H. Kodama, A. Ishibashi, and O. Seto, Phys. Rev. D 62, 064022 (2000). 

[9] S. Mukohyama, Phys. Rev. D 62, 084015 (2000). 
[10] S. Mukohyama, Class. Quantum Grav. 17, 4777 (2000). 
[11] S. Mukohyama, Phys. Rev. D 64, 064006 (2001). 
[12] D. Langlois, Phys. Rev. D 62, 126012 (2000). 
[13] D. Langlois, Phys. Rev. Lett. 86, 2212 (2001). 
[14] A. Neronov and I. Sachs, Phys. Lett. B 513, 173 (2001). 
[15] C. van de Bruck and M. Dorca, hep-th/0012073, 

[16] H. Bridgman, K. Malik, and D. Wands, Phys. Rev. D 65, 043502 (2002). 

[17] K. Koyama and J. Soda, Phys. Rev. D 62, 123502 (2000). 

[18] K. Koyama and J. Soda, Phys. Rev. D 65, 023514 (2002). 

[19] V. Bozza, M. Gasperini, and G. Veneziano, Nucl. Phys. B619, 191 (2001). 

[20] O. Seto and H. Kodama, Phys. Rev. D 63, 123506 (2001). 

[21] M. Giovannini, Phys. Rev. D 65, 064008 (2002). 

[22] J.D. Barrow and R. Maartens, Phys.Lett. B 532, 153 (2002). 

[23] A. Riazuelo, F. Vernizzi, D. Steer and R. Durrer, preprint hep-th/0205220 (2002). 
[24] S. Weinberg, The Quantum Theory of Fields, Vol I, Cambridge University Press (1995). 
[25] J.H. Taylor, Rev. Mod. Phys. 66, 711 (1994); 

For an up to date review see CM. Will, Living Rev. Relativity 4, 4 (2001): cited on April 2003, 

http : //wwwlivingreviews . org/Articles/Volume4/ 

2001-4will/ 

[26] P. Kocian, Diploma Thesis, Universite de Lausanne (2003). 

[27] N. Straumann, General Relativity and Relativistic Astrophysics, Springer Monograph in Physics 
(1984). 

[28] T. Damour and G. Esposito-Farese, Class. Quant. Grav. 9, 2093 (1992); Phys. Rev. D 54, 1474 

(1996); and Phys. Rev. D 58 042001 (1998). 
[29] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999). 
[30] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999). 

[31] J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, Phys.Rev. D 64 123522 (2001). 
[32] W. Goldberger and M. Wise, Phys.Lett. B475, 275 (2000). 

[33] T.R. Taylor and G. Veneziano, Phys. Lett. B212, 147 (1988); S. Ferrara, N. Magnoli,T.R. Taylor 

and G. Veneziano, Phys. Lett. B245, 409 (1990). 
[34] W. Israel, Nuovo Cimento 48B, 463 (1966). 

[35] J. Stewart and M. Walker, Proc. Roy. Soc. London A341, 49 (1974). 



